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Abstract. Data from KEK on subthreshold p as well as on and 
production in pC, dC and oC reactions at energies between 3.5 and 12.0 A 
GeV are described for the first time within a unified approach. We use a model 
which considers a nuclear reaction as an incoherent sum over collisions of varying 
numbers of projectile and target nucleons. It samples complete events and thus 
allows for the simultaneous consideration of all final particles including the decay 
products of the nuclear residues. The enormous enhancement of the p cross 
section, as well as the moderate increase of meson production in dC and ctC 
compared to pC reactions, is well reproduced. In our approach, the observed 
enhancement near the production threshold is mainly due to the contributions 
from the interactions of few-nucleon groups. 



PACS numbers: 24.10.-i, 24.10.Lx, 25.40.-h, 25.45.-z, 25.55.-e 

1. Introduction 

Subthreshold particle production in a nuclear reaction is understood as production 
below the energy threshold of the considered process in a free nucleon-nucleon (NN) 
collision. It is thus a nuclear phenomenon which may be explained by rather different 
assumptions on the properties of nuclear matter and on the interaction dynamics. 
Hereby, it is an open question to what extent subthreshold particle production is 
governed by properties of the nuclear ground state wave function and to what extent 
by the dynamical properties of nuclear matter, not reflected in the ground state 
description. The problem is far from a final solution at present and evidently requires 
a systematic study of high-momentum transfer processes, among them subthreshold 
particle production. 

First measurements of subthreshold p production in proton-nucleus collisions had 
been carried out [1-4] a few decades ago. Then investigations in nucleus-nucleus 
collisions [5-9] and more recent studies of proton-nucleus reactions [10] followed. 
In [11] for the first time light-ion-induced p production has been investigated and 
an enormous enhancement of subthreshold p production in deuteron-nucleus reactions 
compared to proton-nucleus interactions has been found. This fact has been confirmed 
by the final results of the KEK group [12]. They measured d- and a-induced 
reactions in the energy region between 2.5 and 5.0 GeV/nucleon and observed at 
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3.5 GeV/nucleon p cross sections nearly two and three orders of magnitude larger 
than measured in proton- nucleus reactions. 

Most of the descriptions proposed so far are based on transport calculations [13- 
22], thermodynamical considerations [23-26] or multi-particle interactions [27] and 
are devoted to proton and/or heavy-ion induced reactions. By measuring light- 
ion-induced subthreshold p production the KEK group [12] intended to provide an 
experimental test for models in this transition region between proton and heavy-ion- 
induced reactions. To the best of our knowledge only one paper [19] has been published 
so far, which compares proton- and deuteron-induced p production at subthreshold 
energies. 

It should be stressed that all mentioned papers on subthreshold p production 
consider the p spectra without any relation to results concerning other reaction 
channels. At KEK [12], however, the spectra of 7r ± and K ± mesons were measured 
together with those of the antiprotons. In a recent paper [28], in the framework of 
the Rossendorf collision (ROC) model we simultaneously considered all these reaction 
channels for pC reactions and achieved a good reproduction of the data. In this paper, 
we extend these considerations to dC and aC reactions using the same parameter set. 

In the ROC model, it is assumed that the nuclear residue becomes excited during 
the reaction due to the distortion of the nuclear structure by the separation of the 
participants from the spectators and due to the passage of the reaction products 
through the spectator system. In this way, final-state interactions are taken into 
account without making special assumptions concerning re-absorption, re-scattering, 
self-energies, potentials etc for each particle type (see also [28]). The fragmentation of 
the residual nuclei and the interaction of the nuclcons participating in the scattering 
process are treated on the basis of analogous assumptions. Even more important, 
the phase-space of the complete final state consisting of the reaction products of the 
interaction of the participants as well as the fragments of the decay of the spectator 
systems is calculated, exactly This feature seems to be unique to the ROC model. 

The plan of the paper is as follows. In sectional the main ingredients of the ROC 
model are explained, which is used for the calculations to be presented. Section [3] 
contains a comparison of theoretical and experimental results for particle production 
with special emphasis on subthreshold p production. A summary is given in section^] 

2. The model 

The ROC model is implemented as a Monte Carlo generator which samples complete 
events for hadronic as well as nuclear reactions. It was successfully tested for pp 
interactions up to ISR energies in [29,30], while nuclear reactions were considered in 
the papers [28,31-34]. A detailed description of the model is given in [30] for the case 
of pp interactions and in a recent publication [28] for hadron-nucleus interactions. In 
the following the basic mathematics from the hadron-nucleus case [28] is extended to 
nucleus-nucleus reactions. Although in this paper we consider only light-ion-induced 
reactions, the given presentation is applicable to any nucleus-nucleus reaction. For 
the interested reader some minor changes in the ROC model relative to previous 
publications are indicated where appropriate. 

The cross section of the interaction of two nuclei (A,Za) and (B,Zb) 
characterized by their mass and charge numbers is considered as an incoherent sum 
over contributions from varying numbers a and b of nucleons (thereof z a and Zb 
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protons) participating in the interaction 

A min(a,ZA) B minibus) 

Ms) = Y. E E E 

a=l z a =max(0,a— Ma) b=l z b =max(0,b—MB) 

dW(s;a aZabZb ) 

Here, s = P 2 denotes the square of the centre-of-mass energy of the projectile-target 
system with P — (E, P) being the total 4-momentum. The summation extends over 
all possible charge states ranging from clusters consisting completely of neutrons to 
those composed totally of protons. This range is restricted in cases where the number 
of available protons Z A (Zb) or neutrons Ma = A—Za (Mb = B — Zb) is smaller than 
the number of cluster nucleons. Expression a aZa bz b stands for the partial cross section 
of the interaction of clusters (a, z a ) with (b,Zb) , and the quantities dW(s; a a z a bz b ) 
describe the relative probabilities of the various final channels OL az j, Zb . 

The partial cross sections a a z a bz b are calculated using a modified version of a 
Monte Carlo code [35] which is based on a probabilistic interpretation of the Glauber 
theory [36]. They account for sequential collisions between a projectile and b target 
nucleons capable of sharing their energy in close analogy with the virtual clusters of 
the cooperative model [37-44] (see [28] for a discussion of the novel features of the 
ROC model compared to the cooperative model). We use the profile function 



T A B(d) 



A B 



i-nn^-^) 



i=ij=i 

A B 



n^(**<) d3r «n^ B ^) d3r i 

of the considered nuclei, which depends on the nucleon densities p A (^i)i Psi^j) and 
the probability 

Pij = exp(-d| 7 -7r/(7jvAr) 

for an interaction of the ith projectile and the jth target nucleon with dy being the 
distance between the interacting particles. The nucleon density [45] 

PA (r) ex (1 + 77[1.5(/ 2 - e 2 )// 2 + e 2 r 2 // 4 ]) exp(-r 2 // 2 ) (2) 

of light nuclei 2 < A < 20 can be derived from a standard shell model wavcfunction 
with rj = (A - 4)/6 , e 2 = / 2 /(l - I /A) and / = 1.55 fm. For the deuteron the Paris 
deuteron wavefunction [46] is used. Then the NN cross section (Jnn is adapted such 
that the integral of the profile function over the impact parameter d reproduces the 
total inelastic AB cross section 



o-% = J d 2 dr AB (d). 



The same calculation yields also the partial cross sections o a z a bz h we are interested in 
(for further details see [35]). 
In the quantities 

dW(s; a.az a bz b ) oc dL n (s; ac aZabZb )p A (P A )p B (PB)T 2 (3) 
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describe the relative probabilities of the various final channels OL aZabZb . They are given 
by the Lorentz-invariant phase-space factor dL„(s; ot a z a bz b ) multiplied by the square 
of the empirical reaction matrix element T 2 describing the collision dynamics. The 
Fermi motion is implemented via the momentum distributions of the residual nuclei 
Pa(Pa) and pb(Pb), which are functions of the numbers a and b of participants. They 
are taken as Gaussians having, in case of nucleus A, a width of 



a a = y/aA/5/(A-l) Pf (4) 

in accordance with the independent particle model [47] with A = A — a being the 
number of nucleons in the residue and pp is the Fermi limit of the nucleus considered. 
For nucleus B an analogous formula holds, and for the deuteron again the Paris 
deuteron wavefunction [46] is employed, now in momentum space. No special high- 
momentum component [48-53] is used in this paper. 

The Lorentz-invariant phase-space is defined as the integral over the 3-momenta 
of the n primarily produced final particles with energy-momentum conservation taken 
into account 

dL n ( s;aaZabZb ) = TTl£*a4(P -y>), (5) 

2—1 I— 1 

Here, the 4- momentum of the ith particle is denoted by pi — (e,, pi) with pf = m 2 . 

For numerical calculations, the <5-function in equation JSJ has to be removed by 
introducing a new set of 3n — 4 variables to replace the 3n 3-momentum components. 
It is reasonable to choose a set of variables, which reflects the underlying physical 
picture of the interaction process. Using recursion [54] equation can be rewritten 
in the form 

dW{s;a aZabZb ) oc -^^P^ P a)^^Pb{Pb) 

&W a {M 2 a )&Wb{M b )&W c {mI) (6) 

with the 4- momenta of the nuclear residues Pa = {Ea,Pa), Pb = (Eb,Pb) and of 
the participants 

Pc = P-Pa-Pb = (E c ,Pc), (7) 
the corresponding invariant masses are given by M\ = P\, M B = P B and Mq = 
P%. The integrals over the Fermi motion d 3 PaPa{Pa) j1E A and $P b pb{Pb)I^E b 
separate, in accordance with the participant-spectator picture, the phase-spaces of the 
riA and ng nuclear fragments 

dW A (M 2 A ) = &M 2 A&L nA {M 2 A )T 2 A dW B {M 2 ) = dM 2 dL nB (M 2 )T 2 (8) 

from the phase-space 

dW c (M 2 ) = dL nc (M 2 )T 2 (9) 

of the nc = n — tla — n B final particles arising from the participant system. In (JHJ 
and 10, the matrix element T 2 — T\T B T^ is split into factors describing residue 
fragmentation {T\ and T B ) and participant interaction (Tjy). There is, however, 
a strong kinematic link between participants and spectators, since invariant mass 
Mc of the participant system, invariant masses Ma and M B of the residues and 
the relative kinetic energy y/s — Ma — M B — Mc of these three particle groups are 
connected by energy-momentum conservation. The excitation energy of the target 
residue comes usually into play via the spectral function (see e.g. [50]) derived from 
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electron scattering data. Unique to the ROC model, however, is the treatment of the 
spectator system in close analogy with the participant subsystem, since it calculates 
the complete final state of both the participants and the spectators. 

The term dWc(M^) in ijfJJl and © describes the interaction of the groups of 
participating nucleons. Such a cluster-cluster reaction is treated in complete analogy 
to a hadronic reaction. In a first step intermediate particle groups called fireballs 
(FBs) are produced, which decay into so-called primary particles. The primary 
particles define the channels for which the weights Q are calculated. Among them 
are resonances, which decay subsequently into stable hadrons. The dynamical input of 
the cluster-cluster reaction is implemented by the empirical transition matrix element 

Tl = T?T* S T* X T?T?T* , (10) 

which describes the interaction process Tf resulting in the production of N > 2 FBs, 
the production of hadrons T^ s via the creation of quark-anti-quark (qq) pairs , the 
invariant-mass distribution of the FBs T^ x , the transverse Tj 2 and longitudinal T 2 
momentum distribution of the FBs, and finally, some factors T 2 t necessary for the 
calculation of the statistical weights. The interaction is assumed to proceed via colour 
exchange leading to the removal of valence quarks or of gluons from the interacting 
hadrons. Additional up, down and strange quark pairs are created in the ratio 

u : d : s = 1 : 1 : X s (11) 

with A s = 0.15. They form the varying number of FBs, which subsequently decay 
into the final hadrons. The transverse momenta Pt of the FBs are restricted by an 
exponential damping (longitudinal phase-space) according to 

N 

T t 2 = []exp(~ 7 P u ) (12) 
i=i 

with the mean p t = 2/7. Two leading FBs, the remnants of the incoming clusters, 
get in the mean larger longitudinal momenta than the central FBs by weighting the 
events with 

Ti = {X x X 2 f. (13) 

Here, the light-cone variables X\ — (E% + P z , x )/(e a +p Zia ) and X 2 — (E 2 — P Zt i)/{zb — 
p z> b) are used with the 4-momenta of the participants given by p a — (e a ,p a ) and 
Pb = (e&,p*b). In [30], the high-energy approximations 

X l = (Ei +P z ,i)A/£ and X 2 = (E 2 - P.^/y/i 

were applied for the more symmetric pp interactions without any substantial influence 
on the results. For completeness, it should be mentioned that in the older publication 
[29] a transition matrix element with transverse damping only for the central FBs 

N 

T t 2 = nexp(- 7 (Pt,i) 2 ) 
1=3 

was applied while the 4-momenta t a i and tb 2 transferred to the leading FBs were 
restricted via 

T 2 = exp 03 (t al +t b2 )). 
This variant worked well up to energies of about 50 GeV. 



Subthreshold antiproton production in pC, dC and aC reactions 



G 



Each FB is characterized by two parameters, a temperature Q FB an d a volume 
Vpb ■ The temperature determines the relative kinetic energy of the particles the FB 
decays into via 

N 

TU®fb) = Y[{M I /e FB )K 1 (M I /Q FB ), (14) 
1=1 

while the volume defines the interaction region and influences mainly the particle 
multiplicity via the statistical factor T^ t oc V FB ~ X ■ In l|14fl K\ stands for the 
modified Bessel function. Final hadrons are built-up by random recombination of the 
available quarks during the decay of the FBs. This procedure automatically ensures 
the conservation of all internal quantum numbers. Subsequently, resonances decay 
until the final state consisting of stable particles is reached. For a more detailed 
discussion of the hadronic matrix element the reader is referred to [30]. In this paper, 
we use the same set of parameters as in [30] for the description of the interaction of 
a single projectile nucleon with a single target nucleon. Most of the other terms in 
JQl, however, contain clusters consisting of several nucleons. The basic parameters of 
the FBs emerging from such an interaction are fixed by scaling the volume and the 
temperature parameter according to 

V FB =V^ B (a + b-l) and Q FB = Q^% x (a + b - l)" 1 / 3 . (15) 

It remains to consider the target residues, the structure of which is strongly 
disturbed by the interaction of the participants and subsequent final-state interactions. 
This leads to the excitation and decay of the spectator systems, which are characterized 
by the two parameters temperature Q B and volume Vr (R — A or B in dependence 
on the residue under consideration) in the same way as the FBs emerging from the 
participant interaction. The part of the matrix element responsible for the residue 
fragmentation 

Ifl(a&.taJ = Tl(<3 R )Ti (a* abZb ) R — A,B 
is identical with the corresponding factors l|10(l applied to the hadronic FBs. In order 
to restrict the excitation energy transferred to the residue we use the asymptotic 
approximation of (I14fl for large mass and small temperature 

Tl(Q R ) = y/M R /Q R exp(-M R /e R ) R = A,B. (16) 
An impact parameter dependence is assumed for the temperature parameter. 
This seems to be reasonable, because a peripheral collision with only few participating 
nucleons should excite the nuclear residue much less than a central collision with many 
participants. As a first guess we use 

A - e«»*[l _ eocp (-a/oA 1 / 8 ) e B = e i S ax [l-cxp(-b/aB 1 ^) (17) 

with a as parameter, here fixed to a — 0.5, which determines how fast the maximal 
temperatures 0™ ax and Q T ^ ax are reached with increasing number of participants. 
The value 0™ ax = Q n B lax = 12MeV has been fixed in [28] from a comparison of 
the energy spectra of fragments from the reaction of 2.1 GeV protons with carbon 
measured at 90° [55] with ROC-model calculations. 

All factors still necessary for a correct calculation of the relative weights of the 
various channels are collected in complete analogy with the corresponding factors for 
the FBs in the product of two terms with R = A and R = B according to 

/ y \ n R 

T£(<xf Zab J = g(a« ab J ^ H&ai + 1)2™,. (18) 
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It contains the spin degeneracy factors (2<7j + 1) and the volume Vr in which the 
particles are produced. The quantity g(ot^ Zab ) is the degeneracy factor for groups of 
identical particles in the final state of the residue decay and prevents multiple counting 
of identical states. Further details of residue treatment can be found in [28]. 

3. Comparison with experimental data 

In figure Q the momentum spectra of 7T , and p from pC, dC and aC reactions 
measured at KEK [12] are displayed. The experimental efforts of the KEK group were 
motivated by measurements of large cross sections of subthreshold p production in 
nucleus-nucleus reactions at LBL-BEVALAC [6,7] and at GSI [8]. The authors [12] 
claimed that the use of light-ion beams for investigating subthreshold p production 
should be useful to verify theoretical models. Secondary effects like re-scattering, 
re-absorption, self-energies, potentials etc., which hide and influence features of the 
primary production process, should be smaller than in heavy-ion reactions. Indeed, it 
is the outstanding feature of the data that in d- and a-induced reactions an enormous 
enhancement of p production by nearly two and three orders of magnitude compared 
to p-induced interactions could be observed. The very fact of finding this enhancement 
in light-ion induced reactions makes this interaction type a candidate for the key to a 
deeper understanding of subthreshold p production in nucleus-nucleus reactions. 

The KEK group [12] interprets their p data by using the 'first-chance NN 
collision model' from [48] where the internal nuclcon momenta were extracted from 
backward proton production [49] as a superposition of two Gaussian distributions. 
In this way, the momentum spectra and the incidence-energy dependence of p- 
induced reactions could be successfully reproduced by adapting one normalization 
parameter. However, the application of this model to d-induced reactions using 
different deuteron wavefunctions and the normalization from p-induced interactions 
severely underestimates the p cross sections at subthreshold energies. Thus, according 
to [12], the effect cannot be explained by the internal motion of the nucleons in the 
deuteron, even if a deuteron wavefunction with a high-momentum component is used. 

In [19], p production in p+A and d+A reactions is analysed within a phase-space 
model incorporating the self-energies of the baryons. It is claimed that the internal 
momentum distribution of the deuteron provides a natural explanation of the large 
enhancement under discussion. 

To the best of our knowledge there is no paper which describes the further increase 
of subthreshold p production in a-induced reactions. Thus, the approach presented 
here seems to be the first attempt to consider the whole set of projectiles in a unified 
picture. Beyond it, we also regard the whole set of ejectiles measured. This is of 
special importance since there is not only the enhancement of p production, but also 
the increase of the K~ cross sections with increasing energy and mass number of 
the projectile around the elementary production threshold at 2.6 GeV. And also the 
completely different energy and projectile mass dependence of the pion production 
cross sections far above the threshold should be explained by a realistic approach. 

In figure ^ the results of the ROC model calculated with one fixed parameter 
set are compared to the data [12]. The overall agreement is quite satisfactory in view 
of the different projectile types, the large region of incidence energies, the variety of 
ejectile species and the huge differences of many orders of magnitude in the considered 
cross section values. Particle yields are influenced by the suppression factor A s = 0.15 
(equation ifTTjO of strange quarks and by the algorithm for creating the final hadrons 
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momentum p (GeV/c) 



Figure 1. Momentum spectra of 7T , K and p [12] (symbols) from pC, dC and 
aC reactions compared with ROC-model calculations (histograms). Experimental 
and calculated results for 7r+ and K~ mesons are multiplied by the factors 
indicated in the legend. The scale extends from 10 -8 to 10 4 mb sr — 1 GeV - 1 c 
for energies up to 4.0 A GeV, while the lower limit is 10 -6 for 5.0 A GeV and 
10" 2 mbsr- 1 GeV~ 1 c for 12.0 A GeV. 
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from the quarks produced in the first stage of the interaction process. Hadrons are 
built-up in each FB independently according to the rules of quark statistics [56] by 
randomly selecting sequences of quarks q and antiquarks q. A qq gives a meson, while 
baryons or antibaryons are formed from qqq or qqq. From a given sequence of quarks 
the different hadrons are formed according to the tables of the particle data group [57] . 
There is no parameter which directly determines the ratio between meson and baryon 
production as, e.g., in the PYTHIA-LUND model [58-61]. Only an indirect influence 
via the temperature and the volume parameters is possible, which changes the relative 
weights of the FBs in dependence on their invariant mass and final particle multiplicity, 
respectively. (See [30] for a more detailed consideration of hadron formation.) 

It is, therefore, remarkable that the general trend of the energy dependence of 
the data for all projectile types, moderate increase with energy of the cross sections 
for the light mesons and steep increase of that for the antiprotons, is well described 
by the calculations. Also the shift of the maximum in the p spectra towards higher 
momenta is reproduced. A similar tendency can be observed if the cross sections 
from dC and aC reactions are compared to those from pC interactions. The heavier 
the observed particle the steeper the increase of the cross sections for the heavier 
projectiles, especially at the lowest incidence energies. 

In order to make the presented results more obvious we show in figures [21 and 
the partial spectra from collisions of a projectile with b target nucleons 

[a] + lb] 

for the three projectile types and two selected ejectiles p and tt + . These partial spectra 
represent incoherent sums over contributions from all possible charge numbers (sec 
according to 

do-ab(s) E 

z a —max{0,a—J\fA) zj,—max(0,b—AfB) 

dW(s;a aZabZb ) 



Vaz a bzj, 



Ea a2 „ tZl JdW{s;a aZabZb ) ' 



The partial spectra result from the interplay between the values of the integrated 
partial cross sections 



E 

z a — max{Q : a— A/a) z^=max(0,b— A/*b) 



,T »'< : ) E a az a bz b 



and the energy dependence of the relative weights dW(s; a aZa b Zb ) representing the 
production cross section in the [a] + [b] interaction. In the case of p+A reactions the 
crib decrease monotonically with increasing number b of participating target nucleons, 
while for A+A interactions symmetric combinations a — b are preferred. Decisive 
for subthreshold particle production is the steep increase of the production cross 
section with the available energy above the threshold in the [a] + [b] interaction under 
consideration. This explains the small contributions from the [1] + [1] interactions, 
although (Tii represents for all projectile types the largest partial cross section. It is the 
small available energy which prevents higher p production. With increasing number 
of participants both the mean available energy and the number of competing final 
channels increase. This results in maximal contributions from the combinations [1] + [2] 
and [l] + [3] in case of p+A reactions. The spectral distributions are rather independent 
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Figure 2. Partial spectra of p production calculated for pC, dC and aC reactions 
(thin histograms) compared to the data [12] (full circles). In six columns denoted 
by C[b] the spectra from interactions with b participating target nucleons are 
plotted. Likewise, in the rows the numbers of participants in case of d- and a- 
induced reactions is denoted by d[a] and a[a], respectively. The thick histograms 
are the sums of the partial spectra and coincide with the corresponding results 
from figure 111 



of the number of target participants b, see first row in figure [21 This explains the 
successful description of the spectra in the 'first-chance NN collision' approach, which 
corresponds to the [1] + [1] process of the ROC model, by appropriately selecting the 
normalization [12]. Thus, the interpretation of the pC data by the ROC model is quite 
contradictory to the assumptions of a 'first-chance NN collision model'. Comparing the 
[1] + [1] spectra from pC and dC reactions only a moderate increase of p production is 
predicted by the ROC model due to the Fermi motion in the deuteron. The enormous 
enhancement observed stems mainly from processes with both nucleons of the deuteron 
and several target nucleons involved where the energy available for particle production 
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Figure 3. The same as figure 121 for ir + . 



is much larger than in [1] + [1] processes. In case of aC reactions, the number of partial 
spectra contributing essentially increases. Here the larger available energy arises both 
from the larger masses of the interacting subsystems and from their internal motion 
in the projectile and the target. Thus, one can conclude from the viewpoint of the 
ROC model that the key quantity for understanding subthreshold particle production 
is the number of participating nucleons. A direct experimental determination of this 
number for subthreshold p production is highly desirable as discussed in [34] for the 
case of K~ production in pA reactions. 

If we consider particle production far above the threshold as in figure [3] then the 
main difference to subthreshold production consists in the higher value of the energy 
available for the production of the considered particle. This leads to a smoother energy 
dependence of the production cross section for the particle under consideration and 
the gain of available energy due to cluster-cluster processes does not play such an 
important role. In this way, the partial spectra reflect primarily the values of the 
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corresponding partial cross sections a a b with only minor corrections due to the energy 
dependence of the production cross section. Thus, the [1] + [1] processes with the 
largest partial cross section for all projectile types yield also the main contributions 
to the spectra. 

4. Summary 

Subthreshold particle production is a collective phenomenon which is far from being 
completely understood. From the viewpoint of the ROC model, data on subthreshold 
particle production can be reproduced by considering the interaction of few-nucleon 
groups in complete analogy to the interaction of single nucleons, also with regard to 
high-momentum transfer processes. It has been demonstrated here that the cluster 
concept yields a quite natural explanation of the enhancement of subthreshold particle 
production due to the energy gain in the interaction of few-nucleon groups compared 
to NN interactions. This concept should be applicable not only in proton- or light-ion- 
induced reactions, but also for heavy-ion interactions, although in the latter case the 
number of partial processes increases tremendously. In this sense, the ROC model can 
be considered as a promising approach to a unified description of particle production 
processes in a large variety of different types of nuclear reactions. 
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